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The idealization of monochromatic plane waves leads to considerable simplifications in the analysis 
of electromagnetic systems. However, for active systems this idealization may be dangerous due to 
the presence of growing waves. Here we consider a gainy slab, and use a realistic incident beam, 
which is both causal and has finite width. This clarifies some apparent paradoxes arising from earlier 
analyses of this setup. In general it turns out to be necessary to involve complex frequencies uj and/or 
complex transversal wavenumbers kx. Simultaneously real uj and kx cannot describe amplified waves 
in a slab which is infinite in the transversal direction. We also show that the only possibility to have 
an absolute instability for a finite width beam, is if a normally incident plane wave would experience 
an instability. 


I. INTRODUCTION 

Physicists love idealizations and simplifications, to un¬ 
derstand the main principles of physical systems without 
being disturbed by unnecessary complications or degrees 
of freedom. The idealizations may make it possible to 
get analytical results for complicated systems which oth¬ 
erwise call for advanced numerical methods. Without the 
idealizations, the physicist’s understanding and intuition 
would suffer, and the ability to master a large parameter 
space would be very limited. 

However, sometimes the involved idealizations may be 
conflicting or give unphysical results. When consider¬ 
ing the optical or electromagnetic response of a slab, it is 
convenient to let the slab be infinite in the transversal di¬ 
rection, to eliminate complications from reflection at the 
transversal boundaries. At the same time, it is custom¬ 
ary to let the incident wave be monochromatic and plane, 
extending infinitely in time and transversal space. These 
simplifications have been proved very useful for passive 
media, and have been used in a number of settings (see 
e.g. [IH3]). 

Even for gainy slabs, it is quite common to do a 
monochromatic plane wave analysis HE]. Clearly, such 
an analysis is appealing due to its simplicity, and may 
give correct results for certain situations. However, it has 
been argued that this analysis is dangerous in general, 
and may lead to unphysical results Real physics 

happens in the time- and spatial domain, and to justify 
the monochromatic plane wave analysis, Fourier trans¬ 
forms of the fields with respect to time t and the transver¬ 
sal coordinate x must exist. 

In this paper we consider a slab which extends in¬ 
finitely in the transversal x-direction; however we vary 
the description of the source: In Sec. [H] the incident 
beam is described by monochromatic, plane waves. This 
gives peculiar results which we later argue are incorrect. 
In particular, the refraction of the beam inside the slab 
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changes drastically (from positive to negative), subject 
to a small increase of the incident angle. This happens 
despite the fact that the slab is made of a conventional, 
dielectric, weak gain medium. 

It is therefore natural to consider a more realistic 
source. In Sec. m we consider a causal, plane wave 
source, starting at some time t = 0. Due to the possi¬ 
bility of instabilities, we use the inverse Laplace trans¬ 
form to express the time-domain fields from frequency 
components. This analysis, which has been discussed in 
the literature previously EHS], gives physically reasonable 
results. However, it has the disadvantage of predicting 
absolute instabilities related to the infinite extent of the 
plane wave in the x-direction. These absolute instabilities 
are nonexistent in real experiments where the transversal 
dimension of the source is limited m- 

Thus we finally describe the source realistically as a 
causal beam of finite width (Sec. [Tv| ). This analysis, 
although somewhat more complicated, gives the intu¬ 
itive results that we expect from previous experiments. 
Moreover, by deforming the paths of the inverse inte¬ 
gral transforms, it turns out that the system can often 
be described monochromatically, provided we allow the 
transversal wavenumbers k^ to be complex in the super¬ 
position. We derive an explicit expression for the differ¬ 
ence between the solution from the monochromatic plane 
wave analysis and the proper analysis, and demonstrate 
that it is the neglect of this term that lead to the in¬ 
correct results in Sec. [Hj This fact is also independently 
verified by finite-difference-time-domain (FDTD) simula¬ 
tions (Sec. [V[ ). 

Our results are quite intuitive: A gain medium leads to 
growing fields. When a finite width beam has oblique in¬ 
cidence, we expect the fields in the slab to grow along the 
transversal x direction. Therefore, to describe the fields 
in (cj,/ca,)-space, it is necessary to let either the angular 
frequency uj or the transversal wavenumber kx be com¬ 
plex. Modes with simultaneously real uj and kx cannot 
describe amplification in the x-direction in an unbounded 
slab. 

When considering growing waves, an instability is clas¬ 
sified as absolute if the fields grow unlimited with time 
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FIG. 1. The setup consists of a slab with permittivity 62 , 
surrounded by a passive medium ei. The dimensions of the 
slab are inhnite in the x and y directions, with thickness d 
in the z direction. The incident wave is polarized along the 
y-direction. 


even at a fixed point in space mm- Otherwise, an 
instability is convective. The movement of singularities 
of the field as we deform the complex integration paths 
in uj and space [HHH], is used to classify the possible 
instabilities in the slab. It turns out that the only pos¬ 
sibility of having an absolute instability, is if the mode 
kx = 0 is unstable. That is, to find out whether the slab 
supports absolute instabilities for finite width beams, we 
only need to check plane waves at normal incidence, and 
determine if multiple reflections inside the slab are being 
amplified enough to compensate for the transmission out 
of the boundaries. 

The setup is a slab of permittivity 62 surrounded by 
a passive medium ei. Both media are linear, time- 
invariant, isotropic, homogeneous, and without spatial 
dispersion. For simplicity they are assumed nonmag¬ 
netic. The slab is infinite in the x and y directions, and 
has thickness d in the z direction, see Fig. The in¬ 
cident wave is polarized with its electric field along the 
^-direction. 

The linearity assumption needs some clarification. In 
practice any gain medium is only linear for sufficiently 
small fields, due to gain saturation. Thus the linear anal¬ 
ysis is accurate in the absence of instabilities, when the 
incident beam is sufficiently weak. It is also useful to de¬ 
termine the behavior of the system before the instabilities 
have become too large. 


II. MONOCHROMATIC PLANE WAVE 
ANALYSIS 

In this section we assume monochromatic fields of the 
form where uj is real. To compute the electric 

field, we apply a straightforward plane wave expansion 


method, considering the incident beam to be gaussian. 

By solving Maxwell’s equations in the frequency do¬ 
main, and applying the boundary conditions for the elec¬ 
tric and magnetic fields, the electric field amplitude of 
each plane wave component kx is found to be 


E{kx,z,oj) 


'^iKz _^_Re-ik,z^ 2<0 

< 0<z<d (1) 


The total reflection coefficient R, the amplitudes in the 
slab and the transmission coefficient T are functions 
of {kx^uj) through 


{kz -kf)[l - exp{2ik'^d)\ 

(kz + k'^y — {kz — k'y‘^ exp(2i/c^d) 
_ ‘^kzjkz + k'z) _ 

{kz + k'yj‘^ — {kz — k'yy exp(2i/c^(i) 

_ - K) _ 

{kz — k'y)‘^ — {kz + k'y‘^ ex.p{—2ik'^d) 
4:kzk'^ ex.p{ik'^d) 

{kz + k'y‘^ — {kz — k'y‘^ ex.p{2ik'^d)' 


( 2 a) 

( 2 b) 

( 2 c) 

( 2 d) 


Here kz = exuj"^ jc? — k‘^ and k'^ = yje^uj"^ jc? — k‘^ are 
the longitudinal wave number outside and inside the slab, 
respectively, and c is the vacuum light velocity. It is 
worth noting that the expressions 0 do not depend on 
which sign we choose for k '^. The physical electric field is 
given as a superposition of all plane wave components kx^ 
by an inverse Fourier transform. The weighting of each 
component kx is given by the distribution of the incident 
beam. For a gaussian beam the plane wave spectrum is 
gaussian, centered around the main plane wave compo¬ 
nent. 

For concreteness, we consider a slab with permittivity 
given as an inverted Lorentzian, 


62(0;) = 1 - 


Fool 

ujI — uj‘^ — iVuj ’ 


(3) 


where ujq is the resonance frequency, F > 0 (gain) and 
the gain is weak {Fujo/V <C 1). In all calculations here 
and below, we let F = 0.001 and F = O.lcJo. Moreover, 
61 = 2.25, and the thickness d of the slab is given by 
ujodjc = 100. For the monochromatic calculations in 
this section, the observation frequency is set equal to the 
resonance frequency ujq^ at which 62 = 1 — O.Oli. The 
width (7 of the source is given by uj^cTjc = 35. 

Figures and b show the resulting field, for incident 
angles ^ = 30 and 0 = 37 degrees, respectively. A rather 
peculiar observation is made: As the incident angle is 
altered from 30 to 37 degrees, the beam goes from being 
positively to negatively refracted inside the slab. Here 
“negative refraction” refers to the direction of the beam. 
As the slab consists of a weak gain medium with low 
dispersion, this result does not seem to be realistic. 

To investigate this peculiarity further, calculations 
were performed for incident angles from 0 to 90 degrees. 
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FIG. 2. A gaussian beam incident to a gainy slab at incident 
angles of a) 30 degrees and b) 37 degrees. A small change of 
the incident angle cause the refraction inside the slab to go 
from positive to negative. The arrows indicate the direction 
of the phase velocity. 


At small incident angles the transmitted electric field 
is found to increase with increasing 0. Intuitively, this 
makes sense, as the propagation length inside the slab 
increases. However, for incident angles larger than 30 
degrees the field does not increase any more, but instead 
gradually goes over to be negatively refracted. This neg¬ 
atively refracted beam decreases in magnitude with dis¬ 
tance from the entrance of the incident beam. 

As will become apparent from the next sections, these 
results are incorrect, and arise from the assumption that 
the fields are monochromatic and built up of real plane 
wave components. This excludes the possibility of de¬ 
scribing a wave that grows exponentially as x ^ Too. 


III. A CAUSAL SOURCE 

In an actual experiment the source has to be turned 
on at some time t = 0. Due to this abrupt onset, all 
frequencies will be excited to some extent. To determine 
the correct response of the slab we may therefore consider 
the physical electric field E{x, z, t) as a superposition of 
frequency components . Since the slab is active, the 
system may contain instabilities, resulting from the mul¬ 
tiple reflections inside the slab being amplified enough to 
compensate for the transmission out of the boundaries. 
In some cases the electric field is thus not Fourier trans¬ 
formable wrt. t, and interpretation at real frequencies is 
not meaningful. To account for these possible instabil¬ 


ities, the Laplace transform may be used to decompose 
the signal into (possibly) complex frequency components 

[ 6 ]. 

To avoid making the analysis unnecessarily compli¬ 
cated we first consider only a single plane wave compo¬ 
nent kx- Such an analysis of the response of an active 
slab to a causal excitation, where each plane wave com¬ 
ponent is treated separately, has been done in [71 [ 8 ]. For 
concreteness, we consider a source with time dependency 

v{t) = (4) 

i.e. a harmonic oscillation at frequency c^o turned on 
at t = 0 by the Heaviside function H(t). The Laplace 
transform of this function is 


V{u) 



( 5 ) 


which constitutes the frequency distribution of the 
source. The electric field in the time-domain is given 
by the inverse Laplace transform 


1 ni^+OO 

E{x, z,t) = — / dujV{uj)E{k:„ z, 

J 27 — 00 

where E(kx^z^uj) is given by Q. The transform is per¬ 
formed along the line Imuj = 7 , where 7 is some positive 
number larger than the maximal growth rate of the field 
with time. The integration path can be closed along an 
infinite semicircle in the lower half-plane. The residue 
theorem then gives the field as 


E{x, z, t) = E{k^, z, + rem( 2 ). (7) 


Here the first term originates from the pole at cjo, and 
the term rem(t) is a result of all other singularities (cuts, 
poles etc.) of the integrand inside the closed integration 
contour. If all singularities of E{kx^z^uj) are located in 
the lower half-plane Imcj < 0 (for the given kx) these 
transients will decay with time. Then, for sufficiently 
large times, the field is given by the first term in 0 . 

If, on the other hand, the field has poles in the up¬ 
per half-plane, we may still close the integration contour 
along a semicircle in the lower half-plane. However, now 
the term rem(t) will contain terms with frequencies cor¬ 
responding to the poles in the upper half-plane. These 
will grow with time, i.e. we have an absolute instability. 

In the calculations in the previous section, E{kx^z^uj) 
has poles in the upper half-plane for kx corresponding to 
incidence angles larger than Oc = 25.3 degrees. We may 
therefore not assume the field to be monochromatic for 
plane waves incident at ^ assuming real kx in the 

slab. This explains the peculiar refraction in Fig. 

A perfectly plane wave is present everywhere^ and thus 
enters the slab at infinitely large \x\. Plane wave compo¬ 
nents experiencing net round trip gain as they propagate 
inside the slab, will therefore grow unlimited with time 
at any fixed point. In fact, even for a thin slab made of 
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a weak gain medium, these absolute instabilities may oc¬ 
cur for sufficiently large 0. The instabilities are not very 
physical, as they are caused by the assumption that the 
incident wave has infinite width. 

Experiments on amplifying waveguides do not suffer 
from instabilities, even for wave propagation at large 0 
(or small angles with the waveguide axis) [10]. This is 
because in an experimental setup the waveguide has finite 
length and the source does not extend infinitely along the 
waveguide axis. Absolute instabilities due to propagation 
over infinite distances inside the waveguide will therefore 
not be present. To determine whether an unbounded 
slab, such as that in Fig. supports absolute instabilities 
when the source is causal and has finite width, we must 
consider a double Laplace-Fourier transform wrt. uj and 
kx- This is done in the next section. 

IV. A CAUSAL SOURCE OF FINITE WIDTH 

For a causal source of finite width, we find that the 
physical electric field is given by 


ii. As we move uj along the arrow in the upper part of 
Fig. |3}3, the poles will e.g. move as indicated by arrows in 
the lower part. For frequencies along the arrow, we now 
deform the integration path in the complex kx domain 
to avoid these pole trajectories. Now, for kxS along this 
deformed path, there will be no poles for cj’s along the 
arrow, and we may move the encircled segment of the 
cj-integration path down to the real cj-axis. 

To be able to move the encircled uo in the upper figure 
of[^ we required that the integration path wrt. kx could 
be deformed around the pole trajectories. We will later 
show that this is equivalent to requiring that there are 
no absolute instabilities present. If this is the case, we 
may reduce Im uu from 7 to 0 for each segment of the 
inverse Laplace contour. Note that the deformation in 
the complex kx plane in general is different for each uj- 
segment. Equation (§ can thus be rewritten to 

E{x,z,t) = —J dwFMe—*• (9) 

(^tt) J —00 

[ dk^U{k^)E{k^,z,L0)e^'^-^ 

J k{uj) 


2 ni ^+00 

E{x,z,t) = / dwy(w) 

) Jij—oo 


—iujt 


(8) 


f 

J —C 


^kxU{kx^^i.kx't 


ikxX 


by superposition of plane waves of the form Here 
E{kx,z,uj) is given by and U{kx) and V{uj) are the 
plane wave and frequency spectrum of the source, respec¬ 
tively. Note that if the source has finite width, its Fourier 
transform U{kx) is an analytic function (entire function). 

It is reasonable to expect that for a sufficiently weak 
gain medium, or sufficiently small slab thickness, there 
will be no absolute instabilities. If this is the case, it is not 
necessary to describe the field using complex frequencies 
with Imcj > 0. There may however be convective insta¬ 
bilities, i.e. amplification of the field in the ±x-directions. 
This suggests that in certain cases interpretation of the 
fields at real frequencies should be meaningful, provided 
complex wavenumbers kx are included in the superposi¬ 
tion. 

Formally, we will now employ the deformation method 
developed in Ref. [T 2 |, in which the cj-integration in ^ 
is moved down to the real axis, at the expense of de¬ 
forming the /Ca,-integration into the complex plane (Fig. 
[^. To this end we need to know the singularities of the 
field E{kx^ z^uj), which can be found by considering the 
denominator in 0 . In Sec. |III| we noted that even for 
a weak gain medium, the field has poles for Imcj > 0 
for some real kx. These poles suggest that the Laplace 
contour cannot straightforwardly be moved down to the 
real axis. 

However, these singularities may often be avoided, by 
deforming the integration path in the complex kx domain. 
For a given cj, the field will have poles in the complex kx 
plane, for example as indicated in the lower part of Fig. 


where k{uj) is the deformed integration contour in the 
complex kx plane for each uo. The paths k{uo) can be 
written as a sum of integration along the real /c^^-axis, 
plus detours around the pole trajectories. The detour 
integrals may be simplified into contour integrals around 
the poles (lower part of Fig. [^), by noting that the inte¬ 
grations ”up and down” cancel. The integrals around the 
poles are calculated using the residue theorem. For some 
cj’s the poles never cross the real axis, and no deforma¬ 
tion is needed. For other cj’s, several poles may cross the 
axis. In general we get 

E{x,z,t)=— (10) 

J — 00 

(— / dk,U{k,)E{k,,z,uj)0’^^^ 

+ i y] ± Res U{k^)E{k^,z,uj)0^‘^^), 

/c ^ ^ P J 

poles j 


where kpj is the jth pole which has crossed the real kx- 
axis. The sign in front of each residue term is positive 
for anti-clockwise and negative for clockwise integration 
contours. 

We have thus found that in the absence of absolute 
instabilities, ^ may be rewritten into (10). For the 
excitation may be evaluated by closing the uo- 

integration contour along a semicircle in the lower half¬ 
plane. For sufficiently large times, when all transients 
have died out, the electric field is given by 


F(x, 2 ), t) 


1 

27r 



dk,Uik,)E{k,,z,ujo)0’^^^-^^°* 


i y] ± Res 

^cc —^PJ 

poles j 

( 11 ) 
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FIG. 3. (a) The original integration paths in Eq. §. In the lower figure the X’s indicate poles of E(kx, z, uj), for the uj encircled 
in the upper figure, (b) The poles will move as indicated by arrows in the lower figure, as we reduce Imcj from 7 to 0 for the 
encircled uj. We deform the kx integration path for this u) to allow this movement, (c) Provided such a deformation is possible 
for all u) along the integration path in the upper figure of (a), the inverse Laplace transform can be taken along the real axis. 
The lower figure shows the deformed inverse Fourier transform integration path wrt. kx in the complex kx domain, for the 
encircled u in the upper figure. 


In (11) the first term describes the method in Sec. [Ilj so 
the difference is the sum of residues. This sum goes over 
poles that have crossed the real /c^^-axis as Imcj was re¬ 
duced down to 0. These complex plane wave components 
will cause an exponential growth in the ±x-directions, 
and thereby describe convective instabilities. 

The pole trajectories may be found numerically, by 
considering \T{kx,uj)\ in the complex k^ plane for 00 along 
the arrow in the upper part of Fig. Ilf- For the setup 
in Sec. |n[ a total of 34 (17 for Re/c^, > 0 and 17 for 
Kekx < ^ poles cross the real kx-dixis as we reduce Imcj 
from 7 to 0 for Recj = cjo- Figure]^ shows the resulting 
field when the sum over residues has been added to the 
field from Fig. 

The beam propagation in Figs. Ef and If are fun¬ 
damentally different. The summation over residues has 
canceled out the seemingly negative refracted beam in 
If, and the field instead propagates in the +x-direction, 
growing in amplitude. The refraction pattern in Figs. 
If and if are virtually identical, so the summation over 
residues did not affect the field significantly in this case. 
This happens despite the fact that the same observation 
frequency is used, so the summation is over the same 
residues in both cases. The explanation for this is that 
the exponentially growing residue terms are excited dif¬ 
ferently, due to different angles of incidence. The poles 
are located in the region of Re/c^^’s corresponding to an¬ 
gles of incidence between 34.3 and 41.7 degrees. 

In the deformation of the integration paths in ^ to 
the ones in ( 10 ), we required that it is possible to deform 


the inverse Fourier transform path wrt. kx^ to avoid the 
poles crossing the real /c^^-axis as we reduce Imcj from 7 
towards 0. Such a deformation is impossible if two poles 
in the complex -plane collide in such a way that the 
integration path gets “stuck” (see Fig. D). It is argued 
in Appendix that kx = 0 is the only plane wave com¬ 
ponent at which this can happen. If two poles collide 
8it kx = 0 for some Imcj > 0 , we cannot move this uj 
further down, and the inverse Laplace transform has to 
contain complex frequencies with Imuj > 0. The field will 
therefore grow with time, even at a fixed point (absolute 
instability). 

It is perhaps not entirely obvious that the Fourier 
transform in x exists, provided the Laplace transform 
in t is used. To justify ([^, we can use one of the follow¬ 
ing two approaches: Either we just assume the existence 
of the Fourier-Laplace transform, and verify that the so¬ 
lution in the time-spatial domain is consistent with the 
assumptions. Or we consider the transforms in the oppo¬ 
site order. By causality, E{x^ 2 :, t) has finite support as 
a function of x, for any fixed t. Thus it is Fourier trans¬ 
formable. Interchanging the order of integration can be 
done in light of Fubini’s theorem[T5], which places con¬ 
ditions on the source U{kx)V{uj). 


V. FDTD SIMULATIONS 

The analysis and calculations so far have been done in 
the frequency domain, as Maxwell’s equations are easy to 
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FIG. 4. The electric field for the same situation as in Fig. 
but where the summation over residues (the second term in 
has been added to the field. The arrows indicate the 
direction of the phase velocity. 


FIG. 6. FDTD simulations of the same setup as in Fig. 
The angles of incidence are (a) 30 and (b) 37 degrees. The 
arrows indicate the direction of the phase velocity and the 
beam propagation. 


Im kx 



FIG. 5. If two poles collide at = 0 the integration path gets 
“stuck”, and we cannot move the u integration path further 
down to the real axis. This means that the field necessarily 
contains complex frequency components with Imct; > 0, so we 
have an absolute instability. 


solve there. It is however crucial to remember that actual 
physics happens in spatial-time domain. It is therefore 
natural to perform finite-difference-time-domain (FDTD) 
simulations, as an independent verification of Fig- 

ure shows the resulting field distribution from FDTD 
simulations of the same situation as in Figs. and 
Clearly, the FDTD simulations agree with the method 
in Sec. |iyi This supports the conclusion that the stan¬ 
dard monochromatic plane wave solution is generally not 
applicable when analyzing gain media. 

The simulations are based on an implementation of a 
FDTD method [16] for dispersive Lorentzian media m- 


The dispersion (§ is therefore taken into account. The 
simulation domain was chosen large, so that the fields 
never reached the boundaries of the domain during the 
duration of the numerical experiment. 


VI. DISCUSSION AND CONCLUSION 


We have argued that the standard description of the 
electric field as a superposition of real frequencies uj and 
real wavenumbers kx^ is not very useful when analyzing 
media with gain. For the situation with an infinite slab, 
such a description requires the electric field to be simul¬ 
taneously Fourier transformable wrt. t and x. This is 
generally not the case, due to amplification in the slab. 
Therefore, kx and/or uj must be complex in the integrals. 
In Sec. IV we derived an expression for the electric field 
in a gainy slab, by assuming a causal, incident beam of 
finite width. Provided there are no absolute instabilities 
present, interpretation of the field at a real observation 
frequency c^o is meaningful. From the explicit expression 
of the field ([TT]), we find that a sum of terms with com¬ 


plex wavenumbers kx must be added to the conventional 
Fourier integral in kx- This fact is also independently 
verified using FDTD simulations. 

From (11) we may further try to approach the plane 
wave limit, by letting cr ^ oo. Any physical source will 
have finite support, meaning the excitation is identical 
to 0 for sufficiently large x. It can be shown that for any 
u{x) with finite support, its Fourier transform U{kx) will 
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diverge for complex kx^ when cr ^ oo. This means that 
the summation over residues in ( 11 ) will lead to diverg¬ 
ing field in this limit. In the monochromatic limit, the 
plane wave limit a ^ oo is thus not meaningful. The 
physical explanation for this is that in the limit a ^ oo, 
the field from the source propagates an infinite distance 
along the x-axis before reaching a given point, thus pick¬ 
ing up an infinite amount of gain. If there are no poles 
to sum over in 0 (as always will be the case for pas¬ 
sive media), the limit cr ^ oo exists, and the standard 
solution E{x^z^t) = E{kx, is obtained. 


For a realistic incident beam, which is both causal and 
has finite width, we prove that absolute instabilities only 
occur if the kx = 0 mode has poles in the upper half¬ 
plane Imcj > 0. In other words, to check if the slab 
supports absolute instabilities, it suffices to consider a 
normally incident plane wave, and determine if it leads 
to a diverging field. When the kx = 0 mode experiences 
an instability, the slab will support absolute instabilities 
even for oblique incident beams of finite widths, as a 
result of their Fourier spectrum containing the kx = 0 
mode. 


Given that the refraction pattern in Fig. is incor¬ 
rect, one may ask if a similar pattern could arise in a 
more complicated medium. From (11) it is seen that the 
field distribution would in fact be as in Fig. if there 
were no poles to sum over. In other words, if we can find 
a medium e(co’), with e(co’o) = ^ 2 (^ 0 )^ but where no poles 
cross the real /c^^-axis as we reduce Imcj from 7 to 0 , then 
the resulting field would be given by Fig. This will 
not be the case for a conventional weak gain medium, 
such as that described by but will be possible for a 
sufficiently large d using an advanced medium of the type 
considered in [6l [71 [18] . 


Appendix A: Absolute instabilities 


To have an absolute instability, poles of the field must 
collide in the complex /c^^-plane, as we try to deform a cer¬ 
tain segment of the inverse Laplace integration contour. 
This means E{kx^z^uj) must have a pole of order n>2 
in the complex /c^^-plane for a frequency with Imcj > 0 . 
Define 


D = {k,+ k'^f - (k, - k'^f exp(2zA:;d), (Al) 

which is the denominator in (§• Now, for E(kx^z^uj) to 
have a pole of order n > 2 , the derivative dD/dkx must 
be zero, at a zero of D. By differentiation 

= 2{kz + K){^z + K) 

-2{kz-k',){kz-Ky^^'^^ (A2) 

-2ik'^d{kz-k'^fe^^^'^^. 


Here, k'^ mean derivatives of k'^ wrt. kx- From the 
expressions for kz and k'^ we find kz = —kxjkz and k'^ = 
—kxjk'^. The derivative is evaluated at a zero [D = 0), 
so we may simplify ( |A 2 [ ) to 

^ = (A3) 

LlAij, rvz'^z 


From (Al) we see that there are no zeros for kz + 
k'^ = 0. The only possibilities for a zero of order n >2 
in the complex /c^^-plane, are therefore kx = d oi kz = 
—2ijd. The last option means that the field is amplified 
in a passive dielectric, which is unphysical for large times, 
when all transients have died out. 
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